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Abstract 

Single-top-quark production via VF-gluon fusion at hadron colliders provides an op- 
portunity to directly probe the charged-current interaction of the top quark. We cal- 
culate the next-to-leading-order corrections to this process at the Fermilab Tevatron, 
the CERN Large Hadron Collider, and DESY HERA. Using a 6-quark distribution 
function to sum collinear logarithms, we show that there are two independent correc- 
tions, of order l/ln(m|/m^) and as- This observation is generic to processes involving 
a perturbatively derived heavy-quark distribution function at an energy scale large 
compared with the heavy-quark mass. 




Figure 1: Single-top-quark production via ly-gluon fusion. 

1 Introduction 

Now that the existence of the top quark is firmly established |JI| , attention turns to testing 
its properties. A powerful probe of the charged-current weak interaction of the top quark 
at hadron colliders is single-top-quark production. The two primary processes are quark- 
antiquark annihilation via a virtual s-channel W boson 0, |^ and ly-gluon fusion, which 
involves a virtual t-channel W boson (Fig. 1) 0, 0. Within the context of the standard 
model, these processes provide a direct measurement of the Cabbibo-Kobayashi-Maskawa 
matrix element Vtb- Beyond the standard model, they are sensitive to new physics associated 
with the charged-current weak interaction of the top quark [0, |, |, 0, [Tl|, |13|, |1^ . 



Both the precise measurement of Vtb and the indirect detection of new physics require an 
accurate calculation of the single-top-quark production cross section. The quark-antiquark- 
annihilation cross section has been calculated at next-to-leading order in QCD, with a the- 



oretical uncertainty of ±Q% [jT5[. The purpose of this article is to calculate the next-to- 
leading-order correction to the H^-gluon-fusion cross section. 

A complete calculation of the next-to-leading order correction to H^-gluon fusion has 



already been presented in the literature |]T6|. However, we show that this calculation is in- 
correct, due to the factorization scheme used to subtract coUinear divergences. We argue that 
the CTEQ 6-quark distribution function used in that calculation although nominally in 



the deep-inelastic scattering (DIS) scheme, is actually not compatible with that scheme, and 
yields incorrect results. To avoid this problem, we perform our calculation entirely in the 
modified minimal subtraction (MS) scheme [|^. Our numerical results differ significantly 



from those of Ref . |T^ . 

We make several other contributions to the calculation of the next-to-leading-order cor- 
rection to the ly-gluon-fusion process: 

1. We show that there are two independent corrections, of order l/ln(m^/m^) and as, 
which are numerically comparable. The leading-order process is qb q't, as shown in 
Fig. 2(a). The 1/ ln(m^/m^) correction is associated with the diagrams in Figs. 2(b), 2(c), 
while the as correction arises from the diagrams in Figs. 3,4. The existence of a correc- 
tion of order 1/ ln(/i^/mQ) is a generic feature of calculations involving perturbatively- 
derived heavy-quark distribution functions at an energy scale /i large compared with 
the heavy- quark mass rriQ. 



1 



(a) 



(b) (c) 



Figure 2: (a) Leading-order process for single-top-quark production, using a b distribution 
function, (b) Correction to the leading-order process from an initial gluon. (c) Subtracting 
the coUinear region from (b), corresponding to a gluon splitting into a bb pair, (b) and (c) 
taken together constitute a correction of order 1/ ln{mf/ml) to the leading-order process in 
(a). 




Figure 3: Order ag correction to the heavy-quark vertex in the leading-order process qb —>■ q't. 
(c) represents the subtraction of the coUinear region from (b). 



2. We perform the calculation in a simple and systematic way using a structure-function 
approach [l^, This allows the calculation to be organized in a straightfoward 



manner, making use of its similarity with deep-inelastic scattering. 

3. We carefully analyze the appropriate scale in the parton distribution functions. We 
show that the correct scale in the light-quark distribution function is /i^ = {Q"^ is 
the virtuality of the W boson), with essentially no scale uncertainty. However, the 
appropriate scale in the b distribution function is yU^ ~ + mf. 

The paper is organized as follows. In Sec. 2 we show that the next-to-leading-order 
corrections are of two types, 1/ ln(mj/m^) and a^. We then argue that these corrections are 
most reliably calculated in the MS factorization scheme. In Sec. 3 we introduce the structure- 
function approach to calculating these corrections. In Sec. 4 we give our numerical results 
and draw conclusions. We give results for the Fermilab Tevatron pp collider for =1.8 
and 2 TeV, the CERN Large Hadron Collider (LHC), a pp collider with = 14 TeV, 
and the DESY ep collider HERA with = 314 GeV. The analytic expressions for the 
next-to-leading-order structure functions are gathered in the Appendix. 
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(a) (b) (c) 



Figure 4: Order «<, correction to the light-quark vertex in the leading-order process qb q't. 

2 Next-to-leading-order corrections 

2.1 1/ ln(m^/m^) correction 

The tree- level diagrams for VT-gluon fusion are shown in Fig. 1. Since the 6-quark mass is 
small compared with rrit, let us neglect it for the moment. If the b quark is massless, the first 
of these diagrams is singular when the final b quark is coUinear with the incoming gluon. 
This kinematic configuration corresponds to the incoming gluon splitting into a real bb pair. 
The propagator of the internal b quark in the diagram is therefore on-shell, and is infinite. 

In reality the b quark is not massless, and its mass regulates the coUinear singularity 
which exists in the massless case. The coUinear singularity manifests itself in the total cross 
section as terms proportional to ln[((5^ + m^)/m^], where = —q^ is the virtuality of 
the W boson of four-momentum q. Since the virtuality of the W boson is controlled by 
the W propagator, is typically less than or of order M^. For readability, we write the 
logarithm as In (m^/m^) in the following discussion (since 3> M^), although we use the 
exact expression in all calculations. 

The total cross section for W^-gluon fusion contains these logarithmically enhanced terms, 
of order asln(m^/mj), as well as terms of order (both terms also carry a factor of a^, 
which we suppress in the following discussion). Furthermore, logarithmically enhanced terms, 
of order a" In" (m^/m^)/n!, appear at every order in the perturbative expansion in the strong 
coupling, due to coUinear emission of gluons from the internal 6-quark propagator. Since the 
logarithm is large, the perturbation series does not converge quickly, and it appears difficult 
to obtain a precise prediction for the total cross section. 

Fortunately, this difficulty can be obviated. A formalism exists to sum the coUinear loga- 
rithms to all orders in perturbation theory |2T|, |2^, ^ . The coefficient of the logarithmically- 
enhanced term is the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) splitting function 
Pqg, which describes the splitting of a gluon into a bb pair. One can sum the logarithms by 
introducing a b distribution function 6(x,/i^) and calculating its evolution with /i (from 
some initial condition) via the DGLAP equations. Thus the b distribution function can be 
regarded as a device to sum the coUinear logarithms. Since it is calculated from the splitting 
of a gluon into a coUinear hb pair, it is intrinsically of order ln(/i^/m^). We elaborate on 
this point at the end of this section. 

Once a b distribution function is introduced, it changes the way one orders perturbation 
theory. The leading-order process is now qb — > q't, shown in Fig. 2(a). This cross section 
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is of order asln(m^/m^), due to the b distribution function (/i ^ nit). The VT-gluon-fusion 
process, shown in Fig. 2(b), contains terms of both order as\n{mf /ml) and a^, as discussed 
above. However, the logarithmically enhanced terms have been summed into the b distribu- 
tion function and thus are already present in Fig. 2(a). It is therefore necessary to remove 
these terms from the IV-gluon fusion process to avoid double counting. This is indicated 
schematically in Fig. 2(c); the double lines crossing the internal 6-quark propagator indi- 
cate that it is on-shell, which corresponds to the kinematic region responsible for the large 



collinear logarithm |21|, |2^, p3| . 

After the subtraction of the terms of order asln(m^/m^) in Fig. 2(b) by the terms in 
Fig. 2(c), the remaining terms are of order a^. Compared with the leading-order process in 
Fig. 2(a), this is suppressed by a factor 1/ ln(m^/mf). Thus the diagrams of Figs. 2(b), 2(c), 
taken together, correspond to a correction to the leading-order cross section [Fig. 2(a)] of 
l/ln(mj/m^), not of order ag- This is an essential point which has been previously over- 
looked. 

This observation is generic to any process involving a perturbatively derived heavy-quark 
distribution function in the region /i^ ^ mq. For example, the calculation analogous to the 
diagrams in Figs. 2(b), 2(c) for charm production in neutral-current deep inelastic scattering 
Il23| corresponds to a correction of order 1/ ln(Q^/m^) for ^ m^. 

Let us elaborate on our contention that the b distribution function is intrinsically of order 
ln(/x^/m^), rather than merely of order a^. If one neglects gluon bremsstrahlung and the 
scale dependence of the gluon distribution function and the strong coupling, one can solve 
the DGLAP equation for the b distribution function analytically [with the initial condition 
6(x,/i2) = at /i = mb] 0, ^ 

where the DGLAP splitting function is given by 

P,,{z) = llz' + {1 - z)'] . (2) 

Equation (|l]) shows that b{x,fi'^) is of order asln{fj? /ml) compared with the gluon distribu- 
tion function. To support this, we show in Fig. 5 the ratio b{x, fi'^)/g{x, ^'^) x 27r/as(/i^) as a 
function of fi for various fixed values of x, using the CTEQ4M parton distribution functions 
p^ . The curves are approximately linear when /i is plotted on a logarithmic scale, indicating 
that b{x,iJ?) (X [as(/i^)/27r] \n{fi^ /ml)g{x, fi"^). 

The b distribution function is on a different footing from the light-quark distribution 
functions. The light-quark distribution functions involve nonperturbative QCD, and must 
be measured (or calculated nonperturbatively). The b distribution function involves energies 
of order mb and larger, so it can be calculated perturbatively; no measurement is necessary. 
Given the gluon and light-quark distributions functions, perturbative QCD makes a definite 
prediction for the b distribution function. 

2.2 as correction 

There are also bona fide ag corrections to the leading-order process qb q't. The diagram in 
Fig. 3(a) is such a correction; it is of order a^ln(m^/m^) [including the factor as\TL{ml/ml) 
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Figure 5: The ratio of the b distribution function to the gluon distribution function, times 
2-71/ a sifi"^), versus the factorization scale fi, for various fixed values of x. The curves are 
approximately linear when fi is plotted on a logarithmic scale, indicating that b{x, iJ?) oc 
[a<j(/i^)/27r] ln(yU^/m^)(yf(a;, /i^), as suggested by the approximation of Eq. (|l|). 

from the h distribution function], so it is suppressed by a factor of as with respect to the 
leading-order process. 

The diagram of Fig. 3(b) contains terms of both order aj. ln^(m^/m^) and a1 ln(m^/m^). 
The former terms arise from the coUinear emission of the gluon, which gives rise to another 
factor of ln(m^/m^) (on top of the factor from the h distribution function). Similar to 
the discussion above, another power of this logarithm appears at every order in the strong 
coupling, and summation is required to improve the convergence of perturbation theory. The 
coefficient of this logarithmically enhanced term is the DGLAP splitting function P^g, which 
describes the splitting of a quark into a quark and a gluon. The collinear logarithms are 
summed by adding another term, corresponding to gluon emission, to the DGLAP evolution 
equation for the h distribution function. Once this is done, the collinear region must be 
subtracted from Fig. 3(b); this is shown schematically in Fig. 3(c). The remaining terms are 
of order a^ln(mf/mf), so they are bona fide corrections to the leading-order process. 

Finally, there are the corrections to the light-quark vertex in the leading-order process, 
as shown in Fig. 4. These are also bona fide corrections. Figs. 4(a), 4(b) contain collinear 
logarithms \n{Q'^ /m?^ (where niq is a light-quark mass) which are absorbed by the light- 
quark distribution functions in the usual way. Since the light-quark distribution functions 
are intrinsically of zeroth order in a^, the remaining corrections are of order a^. 
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Figure 6: Next-to- next-to-leading-order contribution to single-top-quark production via W- 
gluon fusion. 

2.3 Higher orders 

Consider the next-to-next-to- leading order diagram in Fig. 6. This diagram generates terms 
of order ln^(mf/m^), al\n(mf/ml), and a^. The term of order In^ (m^/m^) comes 
from the region in which the initial gluon splits into a collinear bb pair, and the b quark 
subsequently radiates a collinear gluon. This term is summed by the leading-order DGLAP 
equation, which sums leading logarithms {m^ / ml) / n\ , as discussed in Sec. 2.1. Thus 

this term is already present in the leading-order diagram. Fig. 2(a). 

The terms of order a^ln(m^/m^) come from two sources. The first is when the initial 
gluon splits into a collinear bb pair, and the b quark subsequently radiates a noncollinear 
gluon. This is associated with the diagrams in Figs. 3(b) and 3(c), taken together, which 
correspond to noncollinear gluon radiation. The logarithm is summed via the leading-order 
DGLAP equation into the b distribution function in Figs. 3(b), 3(c), so this term is already 
accounted for. 

The other term of order a^ln(mj/m^) is summed by extending the DGLAP splitting 
function Pqg to next-to-leading order. This sums the first subleading logarithms, of or- 
der \n^~^ {mf / ml) {n > 2) into the b distribution function of the leading-order process. 
Fig. 2(a). The remaining term, of order af, is a correction of order x l/ln{ml/ml) 
compared with the leading order process of Fig. 2(a). 

This analysis demonstrates that all collinear logarithms are ultimately summed into the 
b distribution function; no explicit collinear logarithms remain. The remaining terms are all 
of order a" or, if the diagram has a b quark in the initial state, of order a" ln(mj/m^). These 
correspond to corrections of order a""^ x l/ln(m^/mf) or respectively, compared with 
the leading-order process. For a more detailed discussion of higher orders, see Ref. . 



2.4 Factorization scheme for heavy quarks 

The factorization scheme used to eliminate the collinear divergences from the parton cross 
section must be the same as the scheme used to define the parton distribution functions in 
order to yield a correct (and scheme-independent) result. In the MS scheme, the b distribu- 
tion function b{x, /i^) is defined to be zero at /i = mj,, and is then evolved to higher values 
of n via the DGLAP equations This is the definition of the b distribution function 
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Figure 7: Single-top-quark production via ly-gluon fusion from a structure-function point 
of view. The W boson initiates deep inelastic scattering on both hadrons. 



employed in the CTEQ MS parton distribution functions [|17|, |J . 

Another popular factorization scheme is the DIS scheme. In this scheme, the neutral- 
current structure function F2{x, Q^) is defined to have no radiative correction for light quarks. 
For /X ^ mb, the h quark is essentially a light quark, so a natural interpretation of the DIS 
scheme for the h quark is that its contribution to F2{x, Q^) has no radiative correction. This 



is the interpretation that was made in Ref . , which adopted the DIS scheme for the parton 



cross section and used the CTEQ DIS distribution functions |T^. However, the CTEQ DIS 
h distribution function is actually not in the DIS scheme as interpreted in Ref. |]16|. Rather, 
the h distribution function is again defined by the initial condition b{x,fi'^) = at /i = rrib, 
and evolved to higher values of /i via the DGLAP equations. There is no sense in which 
this yields a b distribution function which is formally equivalent to the usual DIS scheme. 
As a consequence, it is not correct to calculate the parton cross section in the usual DIS 
scheme when using the CTEQ DIS b distribution function. The same is true of the CTEQ 
DIS charm distribution function. 

To avoid this problem, we calculate entirely in the MS scheme. This yields very different 
numerical results from the calculation of Ref. ITBi in the DIS scheme. 



3 Structure-function approach 

Inspecting the leading-order process in Fig. 2(a), qb q't, one observes that it is analogous 
to charged- current deep-inelastic scattering. In fact, it is double deep-inelastic scattering; 
the virtual W boson is probing both the hadron containing the b quark, and the hadron 
containing the light quark, q. This is shown schematically in Fig. 7. We can exploit this 
analogy to calculate the corrections to this process in a compact way, in terms of next-to- 



leading-order hadronic structure functions |19, 2C]. This factorization of the process is exact 



at next-to-leading order, because diagrams involving gluon exchange between the light-quark 
and heavy-quark lines do not interfere with the tree diagram, due to color conservation. 

The hadronic tensor describing a W boson of four-momentum q striking a hadron of 
four-momentum P can be written in terms of five structure functions: 

MW,,{x,Q') = F,{x,Q')(-g,, + ^ 
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(3) 



where = —q^. If the struck quark, and the quark into which it is converted, are both 
massless, then the current with which the W interacts is conserved, and one has q^W^y = 
q^W^y = 0. This imphes that the structure functions F4, F5 vanish. The scahng variable x 
is given by x = Q^/2P ■ q, as usual 

If the quark into which the struck quark is converted is massive, such as the top quark, 
then the current is no longer conserved, and F4, F5 are nonvanishing (although we will find 
that they do not enter our calculation). Furthermore, the scaling variable is now given by 
X = {Q"^ + ml)/2P ■ q. 

The hadronic cross section in Fig. 7 is obtained by contracting the hadronic tensors at each 
vertex with the square of the W propagator connecting them. Due to current conservation 
of the light-quark tensor, the q'^q'^/M^^ term in the numerator of the W propagator does not 
contribute, so one simply contracts the two tensors together. One finds 

3Fi{xi,Q^)Fi{x2,Q^) 



+ F2{XUQ')F2{X2,Q'' 



gZ nZ 

1 



Pi ■ qP2 ■ i-q) 
2. / Pi ■ P2q' 




+ Ifs{x,, Q')F,{x2, Q') ( - l] , (4) 



where 



Q' = -q\ (5) 

The heavy-quark structure functions F4, F5 do not contribute to this expression because they 
are the coefficients of tensors which contain g^, g'^, or both. These tensors give vanishing 
contribution when contracted with the light-quark tensor, due to current conservation. The 
W boson interacts with massless quarks in the hadron of four- momentum Pi, and interacts 
with a h quark in the hadron of four- momentum P2, as indicated in Fig. 7. Note that the 
latter hadron is probed hy & W boson of four-momentum — g, which results in P2 ■ (— g) 
appearing in several places in Eq. (^. One must also add the contribution where the W 
boson interacts with massless quarks in the hadron of four-momentum P2 and with the h 
quark in the hadron of four- momentum Pi. 
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The differential hadronic cross section is given by pOU f] 

where VTf = (Pi + g)^ and VTg^ = (P2 — q)"^ are the squared invariant masses of the hadron 
remnants (including the top quark), and S = 2Pi ■ P2 is the square of the hadronic center- 
of-momentum energy. Using 

2Pi-q = W^ + Q\ (9) 
2P2-(-g) = W^ + Q\ (10) 

we can write Eq. (^) in terms of the integration variables Q^, W^, W.^'- 
3Pi(xi,g2)Pi(x2,g') 



where 



The physical region is given by 



Wi > , (14) 
W2 > mt, (15) 
W1 + W2 < 75, (16) 

Q'™^ = ^[5-^^l^-w|±Al/2(5,Wl^w^|)], (17) 

rmn / 

X{a,b,c) = + b'^ + -2ab-2ac-2bc . (18) 

The next-to- leading-order expressions for the structure functions are given in the Ap- 
pendix. We use the MS scheme, for the reasons discussed in the previous section. After the 
subtraction of the coUinear logarithms ln[(Q^ + m'f)/ml], we set the b mass to zero, since 

^This equation is obtained from Eq. (2) of Ref. by setting dT = and integrating out the four- 
dimensional Dirac d function. 
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it is small compared with the top-quark mass.0 When evaluating the next-to-leading-order 
contribution to the cross section, we use the next-to-leading-order expression for the struc- 
ture function corresponding to the light quark or the heavy quark, but not both at the same 
time, as this would yield a contribution of next-to-next-to-leading order. 

Factorization scale 

The similarity of the leading-order process qb — >■ q't with deep inelastic scattering suggests 
that the relevant scale in the light-quark distribution function is /i^ = Q^. If the parton 
distribution functions were extracted solely from deep-inelastic-scattering data at the same 
values of x and relevant to this process, this statement would be exactly correct, because 
the radiative corrections to deep-inelastic scattering are precisely the same as those to the 
light-quark vertex in qb q't. The latter process has additional radiative corrections, both 
to the heavy-quark vertex and between the two quark lines, but these are unrelated to the 
scale in the light-quark distribution function. 

The actual situation is not far from the situation described above. Most of the information 
on the light-quark distribution functions does come from deep inelastic scattering, and the 
relevant values of x and are within the range of the HERA ep collider: x ^ rrtt/ \fS ~ 0.1 
at the Tevatron and x ~ 0.01 at the LHC, with Q'^<My^. We therefore set /i^ = in the 
light-quark distribution function and refrain from varying the scale, as is usually done to 
estimate the theoretical uncertainty from uncalculated higher-order corrections. 

The situation is entirely different for the scale in the b distribution function. The collinear 
logarithm that results from the diagrams in Figs. 2(b) and 3(b) is \n[{Q'^ + ni^t) / ml]. Upon 
subtraction of the collinear region via the diagrams in Figs. 2(c) and 3(c), the remaining 
logarithm is ln[(Q^ + m^)//i^] (see the Appendix). The appropriate scale in the b distri- 
bution function is therefore /z^ ~ + to^. Since the b distribution is obtained from an 
entirely theoretical calculation, we vary this scale in order to estimate the uncertainty from 
uncalculated higher-order corrections. 

The argument above shows that the appropriate scales in the light-quark and 6-quark 
distribution functions are different. Although it may seem unfamiliar to have different scales 
in the parton distribution functions of a given hadronic process, we have shown that it is 
appropriate in this case. The appropriate scale for the production of a quark of mass mq 
via charged- current deep inelastic scattering is = Q"^ + uiq, which yields fj,"^ = for the 
light-quark structure function and /i^ = + mf for the top-quark charged-current structure 
function. 

4 Results and Conclusions 

We evaluate the next-to-leading-order cross section for single-top-quark production via W- 
gluon fusion using the latest CTEQ MS distribution functions, CTEQ4M The cross 

sections at the Tevatron (1.8 and 2 TeV) and the LHC for the sum of t and i production 

^In practice, it is simpler to set the b mass to zero from the outset, and evaluate the cross section in 
N — A — 2e dimensions. The collinear logarithms appear as terms proportional to 1/e — 7 + ln47r, and are 
subtracted in the MS scheme. 
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vs 


LO (pb) 


l/\n{mf/ml) (pb) 


as (pb) 


NLO (pb) 


1.8 TeV pp 


1.84 


-0.39 


0.25 


1.70 


2 TeV pp 


2.67 


-0.55 


0.32 


2.44 


14 TeV pp 


270 


-31 


6 


245 


314 GeV ep 


1.02x10"^ 


-0.34x10"^ 


0.36x10"^ 


1.04x10-^ 



Table 1: Cross sections for single-top-quark production via VT-gluon fusion at the Tevatron, 
LHC, and HERA for rrit = 175 GeV. The cross sections are the sum of t and i production 
at the Tevatron and the LHC, and either t (positron beam) or t (electron beam) at HERA. 
The first column gives the leading-order cross section [Fig. 2(a)]; the second column the 
correction of order l/ln(m^/m^) [Figs. 2(b), 2(c)]; the third column the correction of order 
as (Figs. 3,4); and the last column the next-to- leading-order cross section (the sum of the 
first three columns). All calculations are performed in the MS scheme using CTEQ4M parton 
distributions functions with fi^ = for the light-quark vertex and /i^ = + ttl^ for the 
heavy-quark vertex. 



for^mt = 175 GeV are given in Table 1, assuming Vtb= 1- The leading-order cross sections 
are also evaluated with the CTEQ4M distribution functions. (When evaluated with the 
CTEQ4L leading-order distribution functions, the leading-order cross sections are 1.61, 2.31, 
and 237 pb at the three machines.) The 1 / ln(m^ / ml) and corrections are listed separately. 
The 1/ ln(m^/m^) correction is —20% at the Tevatron, and —11% at the LHC. This confirms 
previous calculations of this correction in the MS scheme [0, ^ |3l|. The as correction is 



+12% at the Tevatron, and +2% at the LHC. The next-to-leading-order cross section is 
the sum of the leading-order cross section and these two corrections. The fact that the 
as correction partially compensates the l/ln(m^/m^) correction is a numerical accident, as 
these are two truly independent parameters. 

Also given in Table 1 is the cross section for e'p Uetb or e~^p Ugtb at HERA 
|p6| , |27| , |28|1 . (The leading-order cross section is 1.21 x 10^^ pb when evaluated with the 
CTEQ4L leading-order distribution functions.) The l/ln(m^/m^) correction is —33%, and 
the as correction is +36%. An integrated luminosity of about 10 fb~^ would be needed to 
produce a single event. This is unattainable given the design luminosity of the machine 
(£ = 1.6 X lO^VcmVs). 

We argued in Sec. 2.3 that the CTEQ DIS b distribution function is incompatible with 
the usual DIS scheme, and yields incorrect results. To demonstrate this, we also perform the 
calculation in the DIS scheme using CTEQ4D distribution functions. The next-to-leading- 
order cross sections at the Tevatron (1.8 and 2 TeV) and the LHC are found to be 2.24, 3.20, 
290 pb. These differ from the results in the MS scheme by much more than the theoretical 
uncertainty in that calculation, which we now estimate. 

To estimate the uncertainty from uncalculated higher-order corrections, we vary the 
scale in the b distribution function about the central value /i^ = + mf. The results 
are shown in Fig. 8 at the Tevatron (2 TeV) and the LHC, for both the leading-order and 
next-to-leading-order cross sections, using the CTEQ4M parton distribution functions. The 
next-to-leading-order cross section is considerably less sensitive to /x, as expected. Varying 



^The current world-average top-quark mass is 175.6 ± 5.5 GeV |2£ 
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H between one-half and twice its central value yields an uncertainty in the next-to-leading- 
order cross section of ±5% at the Tevatron and ±4% at the LHC. As discussed in Sec. 3, we 
do not vary the scale in the light-quark distribution function, where /i^ = Q^. Although our 
estimate of the theoretical uncertainty in the cross section from uncalculated higher orders 
is rather small, it would be worthwhile to pursue the calculation to the next order in a.^- 

Another source of uncertainty stems from the uncertainty in the top-quark mass. The 
cross section as a function of the top-quark mass is shown in Fig. 9 at the Tevatron (2 
TeV) and the LHC. The cross section is relatively insensitive to the top-quark mass because 
the decrease in the parton distribution functions with increasing rrit is not augmented by a 
decrease in the partonic cross section, which scales like l/M^r instead of 1/s. The present 
uncertainty of ±5.5 GeV in the top-quark mass |2^ corresponds to an uncertainty of ±9% 
in the cross section at the Tevatron and ±5% at the LHC. Anticipating an uncertainty of 
±2 GeV in the top-quark mass from Run II at the Tevatron and/or from the LHC reduces 
the uncertainty in the cross section from the top-quark mass to ±3% at the Tevatron and 
±2% at the LHC. 

Another source of uncertainty is the gluon distribution function, which reflects itself in an 
uncertainty in the b distribution function. It is impossible to estimate this uncertainty with 
any confidence at this time; what is needed is a parton distribution set with an associated 
error-correlation matrix. 

In this paper we present the first complete and correct calculation of the next-to-leading- 
order corrections to single-top-quark production via W^-gluon fusion. We show that there are 
two independent corrections, of order 1/ ln(m^/m^) and a^, which are numerically compara- 
ble. We estimate the uncertainty due to uncalculated higher-order corrections to be about 
±5% at the Tevatron and the LHC. Assuming the uncertainty in the gluon distribution func- 
tion can be quantified and reduced to a sufficiently-small level, single-top- quark production 
via VT-gluon fusion will be an accurate probe of the charged-current interaction of the top 
quark at the Tevatron and the LHC. In conjunction with qq —>■ tb, it will yield an accurate 
measurement of Vtb and possibly indicate the presence of new physics. 

Note: The numerical results contained in this paper were obtained by evaluating the 
weak coupling constant g in terms of the Fermi coupling Gp and the ly-boson mass M^y, 
via = 8GfM^/V2, where Gp = 1.16639 x 10"^ GeV^^ and Mw = 80.4 GeV. These 
numerical results are approximately 2% less than the values which appear in the published 
version of this paper [Phys. Rev. D 56, 5919 (1997)]. 
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Appendix 



The structure functions for the charged- current production of a heavy quark were calculated 
at next-to-leading order many years ago in Ref. P^. This calculation was recently repeated 



in Ref. ||33|, which discovered a misprint in the previous result, and also adopted the modern 
convention of treating the gluon as having N — 2 helicity states in N dimensions. We present 
the structure functions below, for completeness. 

Our calculation utilizes the charged-current structure functions for top-quark production, 
Fi{x,Q'^) (i = 1,2,3), calculated in the MS scheme. The bottom-quark mass is neglected 
throughout. To make contact with Refs. ||32|, ^ we define a related set of structure functions, 
J-'i{x,Q'^), via Fi = T\, = IxT^-, and F3 = 2JF3. These structure functions are related to 
the parton distribution functions by 



^f(a;,Q^) = g(a;,/i^)4 



a 



277 



where 



/ - i^f(^,Q^/i^A)g(-,/iM+i^f(^,Q^/i^A) 

J X Z \_ \Z / \Z 



A 



Q'^ + m 



2 ■ 



The coefficient function for real and virtual gluon emission (Fig. 3) is 



r)2 _|_ ™2 

i^f(^,Q^/.^A) = P,,{z)\n^^-^:^ + h'i{zA) 



where0 



h1{z,\) 



4/1 + ^2 



3 V 1 



^\h'^ + Ai5{l-z)+B^,^ 



-B. 



24- 



l-\z) 



3,1 



1-Z 



(1 - xzy 



\ 1 l + 3A^^ \ 



[l + z'^)\nz , 



Ka 



l-z 



-(l-A)ln(l-A) 



21n(l - z)- ln(l - A^) 



l-z 



The coefficients in the expression for h1{z, A) are given in Table 2. 
The coefficient function for initial gluons [Figs. 2(b), 2(c)] is 



Hl,,,{z,Q',^^\X) 



P,a{z) [±L^ + In ^^^) + K{z, A) 



(19) 
(20) 



(21) 



(22) 



(23) 



(24) 
(25) 



(26) 



*The expression for ft,'' corrects a misprint in Ref. [p3|, where the 7r^/3 term was written as tt'^/S. 
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Bo 



B 



3,i 



1 

2 
3 





Ka 




2-2^2-2 

z 



z — z" 
^--1-z 

z 

l-z 



Co A 



3,i 



Table 2: Coefficients in tlie expression for h^z, A). 



1 

2 



wfiere 



4 - 4(1 - A) 

8 -18(1 -A) 
+ 12(1- A)2 

2(1 - A) 



l-Xz 
1-A 



l-Az 



2 

6A 

-2(1-^) 



hl{z,X) 



1 
2 

In 



z^ + (l- zf 



1 - Xz 



1-X)z' 

Co + Ci,iz{l -z)+ C2,i + (1 - X)zLx{C3,i + A^C4,i) , 



Co = Pgp(^) [2 ln(l -z)- ln(l - A^) - In z] 



-4 

-12A 



Table 3: Coefficients in the expression for /if (-2, A). 



(27) 

(28) 
(29) 
(30) 



The coefficients in the expression for /if (z, A) are given in Table 3. 

The explicit logarithms in Hf{z, Q"^, ji^, X) and Hf{z, Q'^, fi^, A) show that the appropriate 
scale for the process is iJ? = Q"^ + ml, as discussed in Sec. 3. 

The structure functions for light quarks (Fig. 4) in the MS scheme can be obtained from 
these expressions by taking — (A — > 1). This limit is unambiguous, except for the 
factor Lx] the correct substitution is 



In 



l-z 



(31) 
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Figure 8: Cross section for single-top-quark production via VF-gluon fusion at the Tevatron 
and the LHC for mt — 175 GeV, versus the ratio of the factorization scale /i to its natural 
value, fjL — ^10^ + m|. Both the leading-order and next-to-leading-order cross sections are 
shown. 
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Figure 9: Next-to-leading-order cross section for single-top-quark production via VF-gluon 
fusion at the Tevatron and the LHC as a function of the top-quark mass. 
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